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' Abstract 

(N 

Q_^| We consider the set multi-cover problem in geometric settings. Given a set of points P and 

(]J • a collection of geometric shapes (or sets) J, we wish to find a minimum cardinality subset of 3^ 

\ such that each point p £ P is covered by (contained in) at least d(p) sets. Here d(p) is an integer 

■ demand (requirement) for p. When the demands d(p) = 1 for all p, this is the standard set cover 

problem. The set cover problem in geometric settings admits an approximation ratio that is 
better than that for the general version. In this paper, we show that similar improvements can 
be obtained for the multi-cover problem as well. In particular, we obtain an O(logopt) approx- 
imation for set systems of bounded VC-dimension, where opt is the cardinality of an optimal 
C/2 \ solution, and an 0(1) approximation for covering points by half-spaces in three dimensions and 

for some other classes of shapes. 
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> ■ 1 Introduction 

The set cover problem is the foUowing. Given a universe U of n elements and a collection of sets 
Q ■ "J = {Si, . . . , Sm} where each Si is a subset of U, find a minimum cardinality sub-collection CCS' 

Q>^ , such that C covers U; in other words, the union of the sets in C is U. In the weighted version each set 

I Si has a non-negative weight Wi and the goal is to find a minimum weight cover C. In this paper, we 

are primarily interested in a generalization of the set cover problem, namely, the set multi-cover 
problem. In this version, each element e G U has an integer demand or requirement d(e) and a 
multi-cover is a sub-collection C C 9" such that for each e G U there are d(e) distinct sets in C that 
^ ' contain e0 The set cover problem and its variants arise directly and indirectly in a wide variety of 

■ settings and have numerous applications. Often 9" is available only implicitly, and could have size m 

exponential in the size of U, or even infinite (for example 3" could be the set of all disks in the plane). 
The set cover problem is NP-Hard and consequently approximation algorithms for it have received 
considerable attention. A simple greedy algorithm, that iteratively adds a set from S' that covers 
the most uncovered elements, is known to give a (1 -|- Inn) approximation, where n = |U|. (In the 
weighted case, the algorithm picks the set with minimum average cost for the uncovered elements.) 
Similar bounds can also be achieved via rounding a linear programming relaxation. The advantage 
of the greedy algorithm is that it is also applicable in settings where 9" is given implicitly, but there 
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exists a polynomial time oracle to (approximately) implement the greedy step in each iteration. It 
is also known that unless P = NP there is no o(logn) approximation for the set cover problem 
|LY94| . Moreover, unless NP C DTIME(n<^(^°g'°s")) there is no (1 - o(l))lnn approximation 
|Fei98j . Thus the approximability of the general set cover problem is essentially resolved if P 7^ NP. 
However, there are many set systems of interest for which the hardness of approximation result 
does not apply. There has been considerable effort to understand the approximability of set cover 
in restricted settings, and previous work has shown that the set cover problem admits improved 
approximation ratios in various geometric cases. In particular, set systems that arise in geometric 
settings are the focus of this paper. 

In the geometric setting, we use (P,9^) to describe a set system (also referred to as a range 
space) where P is a set of points and is a collection of sets (also called objects or ranges). We 
are typically interested in the case where 9" is a set of "well-behaved shapes". Examples of such 
shapes include disks, pseudo-disks, and convex polygons. The goal is to cover a given finite set of 
points P in IR"^ by a collection of objects from 3^. At a higher level of abstraction, one can consider 
set systems of small (or constant) VC dimension. In addition to the inherent theoretical interest in 
geometric set systems, there is also motivation from applications in wireless and sensor networks. In 
these applications the coverage of a wireless or sensor transmitter can be reasonably approximated 
as a disk-like region in the plane. The problem of locating transmitters to optimize various metrics 
of coverage and connectivity is a well-studied topic; see [TWDJOSj for a survey. 

Bronnimann and Goodrich |BG95| . extending the work of Clarkson [Cla93j . used the reweighting 
technique to give an O(loKopt) approximation for the set cover problem when the VC dimension of 
the set system is boundeco- Here opt is the size of an optimum solution. Known hardness results 
[LY94| preclude such an approximation ratio for the general set cover problem. The reweighting 
technique and its application to set cover |Cla93l IBG95| show that the approximation ratio for set 
cover can be related to bounds on e-nets for set systems. Using this observation, [BG95| showed an 
improved 0(1) approximation ratio for the set cover problem in some cases, including the problem 
of covering points by disks in the plane. Long [Lon01| made an explicit connection between the 
integrality gap of the natural LP relaxation for the set cover problem and bounds on the e-nets 
for the set system (see also |ERS05j ). This allows opt in the approximation ratio to be replaced 
by f, where f is the value of an optimum solution to the LP relaxation (i.e., the optimal fractional 
solution). Clarkson and Varadarajan |CV07j developed a framework to obtain useful bounds on the 
e-net size via bounds on the union complexity of a set of geometric shapes. Using this framework 
they gave improved approximations for various set systems/shapes. Recently, Aronov, Ezra and 
Sharir [AES09j . and Varadarajan |Var09j sharpen the bounds of Clarkson and Varadarajan in some 
cases [CyOTl. 

The geometric set cover problem induced by covering points by disks in the plane is strongly NP- 
Hard |FG88j : very recently a PTAS was obtained for this problem |MR09j improving a previously 
known constant factor approximation. Some other geometric coverage problems are known to be 
APX-hard [FMZnTj : that is, there is a constant c > 1 such that unless P = NP, there is no c 
approximation for them. 

Our results. In this paper, we consider the multi-cover problem in the geometric setting. In 
addition to the set system (P,9"), each point p G P has an integer demand d(p). Now a cover 
needs to include, for each point p, d(p) sets that contain p. For general set systems, the greedy 
algorithm and other methods such as randomized rounding, which work for the set cover problem, 

'Bronnimann and Goodrich |BG95| consider the hitting set problem which is the set cover problem in the dual 
range space. In this paper we blur the distinction between set cover and hitting set. 
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can be adapted to the multi-cover problem, yielding a (1 + Inn) approximation (see [VazOlj ). In 
contrast, the e-net based approach for geometric set cover does not generalize to the multi-cover 
setting in a straight-forward fashion. Nevertheless, we are able to use related ideas, in a somewhat 
more sophisticated way, to obtain approximation ratios for the geometric set multi-cover problem 
that essentially match the ratios known for the corresponding set cover problem. In particular, we 
obtain the following bounds. In all the bounds, f < opt is the value of an optimum (fractional) 
solution to the natural LP relaxation, and opt is the value of an optimum (integral) solution. 

• O(logf) approximation for set multi-cover of set systems of bounded VC dimension. 

• 0(1) approximation for (multi) covering points in IR^ by halfspaces. This immediately leads 
to a similar result for multi-cover of disks by points in the plane. 

• O(logloglogf) approximation for covering points by fat triangles (or other fat convex polyg- 
onal shapes of constant descriptive complexity) in the plane. 

The second and third results follow from a general framework for a class of "well-behaved" 
shapes based on the union complexity of the shapes. This is inspired by a similar framework from 
[CV07l[AES09| . Our work differs from previous work for set cover in geometric settings in two ways. 
First, we use the LP relaxation in an explicit fashion in several ways, demonstrating its effectiveness. 
Second, our work points out the usefulness of shallow cuttings for the multi-cover problem. We hope 
that these directions will be further developed in the future. 

2 Preliminaries 

2.1 Problem statement and notation 

Let I = (P,3") be a given set system with VC dimension 6. Here P is a set of points, and 3^ is a 
collection of subsets of P, called ranges or objects. Assume that every point p G P has an associated 
integral demand dj(p) > 0. When the relevant set system is understood, we may write d(p). Here 
we would like to find a minimum cardinality set of ranges of 5" that covers P, such that every p € P 
is covered at least d(p) times. Note that we allow a range of 3' to be included only once in the cover. 
This is an instance of the set multi-cover problem. There is also a weaker version of the problem, 
where the solution may be a multiset; that is, a range may be included multiple times. 

We will also discuss the demand of a set P' C P, which is d(P') = dj(P') = J2peP' "^(P)- The 
total demand of a set system I = (P,9^) is d(P). 

Definition 2.1 For a point p G P and a set X C 3^ where each range in 9" has a non-negative 
weight, let #(p fl X) denote the depth of p in X; namely, it is the total weight of the ranges of X 
covering p. If the ranges do not have weights then we treat them as having weight one. 

Definition 2.2 Given a multiset Z C 3^, let J = (Q,S) = (P,9") \ Z denote the residual set 
system. The residual instance encodes what remains to be covered after we use the coverage provided 
by Z. Each p G P has residual demand drcs(p,-2) = max(d(p) — ^{pCiZ) ,0), and Q comprises 
the points of P with nonzero residual demand. Thus dj{p) = drcs{p,Z). Also Q = 3'\ Z. We will 
also write, for Q' C Q, dres(Q', Z) = X^pgQ/ di.es(p, Z). In particular, dres(Q; Z) = d j{Q) is the total 
residual demand of T, with respect to Q. 

A set system (P, 3^) has VC dimension 6 if no subset of P of cardinality greater than 5 is shattered 
by 5". Here a set P' C P is shattered if for every X C P' there is a range r G 5" such that X = rCiP' . 
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Given a range space S = (P,?"), its dual set system is S* = (9", P*) where P* = {3"p | p € P} and 
IFp = {r G I p € r}. For a set system S with VC dimension 6, we denote by 6* the VC dimension 
of S*. It is known that 6* < 2'^+^ |PA95[ iHarOS] : thus if S has bounded VC dimension, so does S* . 
However, for specific set systems of interest, in particular geometric set systems, one can directly 
show much stronger upper bounds on 6*. 

2.2 LP relaxation 

A standard approach to computing an approximate solution to an NP-hard problem is to solve a 
linear programming relaxation (LP) of the problem and round its fractional solution to an integral 
solution to the original problem. 

In our case, if 9" = {ri, . . . , rm} and P = {pi, . . . , p„}, the natural LP has a variable Xi for range 



m 



min 




i=l 

subject to ^ Xi > d(pj) Vpj G P, (1) 

Xj G [0, 1] i = 1, . . . ,m. 

Note that LP is a relaxation of the integer program for the set multi-cover problem, for which 
Xi are required to take a value in {0, 1}. If repetitions of a set are allowed, then the constraint 
Xi G [0, 1] is replaced by Xi > 0. 

Let f = f(T) denote the value of an optimum solution to the above LP. Clearly, opt > f(I). 
We will refer to the values assigned to the variables Xi for some particular optimal solution to the 
LP as the fractional solution. In the following, we will refer to the value of Xi in the solution as 
the weight of the range r,. We will sometimes use vectors that are not optimal solutions for LP, 
but only feasible; that is, they satisfy the constraints. 

2.3 Overview of Rounding for Geometric Set Cover 

We briefly explain the previous approaches for obtaining approximation algorithms for the set cover 
problem in geometric settings. The work of Clarkson |Cla93| and Bronnimann and Goodrich [BG95j 
used the reweighting technique and e-nets to obtain algorithms that provide approximation bounds 
with respect to the integer optimum solution. In |Lon01l fERSOSj . it was pointed out that these 
results can be reinterpreted as rounding the LP relaxation and hence the approximation bounds can 
also be stated with respect to the fractional optimum solution. Here we discuss this interpretation. 

Note that in the set cover setting d(p) = 1 for all points. Consider a fractional solution to the LP 
given by Xi assigned to ranges G 9", with total value f = "^^Xi. Let e = 1/f. Prom the constraint 
(Eq. ([T])) it follows that for each p, X^i pgr, ^ ^(p) /f = 1/f = e. Interpreting Xi/f as the weight 
of range r^, we obtain a set system in which all points are covered to within a weight of e. Therefore 
an e-net of the (weighted) dual range space is a set cover for the original instance. Now one can 
plug known results on the size of e-nets for set systems to immediately derive an approximation. 
For example, set systems with VC dimension 6 have e-nets of size 0{6/e ■ log 1/e) |PA95j and hence 
one concludes that there is a set cover of size 0((^*f logf) computable in polynomial time, that is, 
an 0{5* logf) approximation. For some set systems improved bounds on the e-net size are known. 
For example, if P is a finite set of points and 5" is a set of disks in the plane then e-nets of size 
0(l/e) are known to exist for the dual set system and hence one obtains an 0(1) approximation 
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for covering points by disks in the plane. Clarkson and Varadarajan [CV07j showed that bounds on 
the size of e-nets can be obtained in the geometric setting from bounds on the union complexity of 
objects in 3^. We remark that the connection to e-nets above also holds in the converse direction: 
for a given set system, the integrality gap of LP can be used to obtain bounds on the e-net size. 

In the multi-cover setting we can take the same approach as above. However, now we have for 
a point p, X^j pgr, ^ d(p) • e where e = 1/f. Note that we now have non-uniformity due to 
different demands and hence an e-net would not yield a feasible multi-cover for the original problem. 

3 Multi-cover in spaces with bounded VC dimension 

In this section, we prove the following theorem. 

Theorem 3.1 Let 2 = (P,3^) be an instance of multi-cover with VC dimension 6 and let 6* be the 
VC dimension of the dual set system. There is a randomized poly-time algorithm that on input T 
outputs 0((5*f logf) sets of 3' that together satisfy T, where f is the value of an optimum fractional 
solution to X. 

We have an easy proof of the above theorem for the setting in which a set is allowed to be used 
multiple times; the proof is based on results on relative approximations. See Section [3?2l for details. 

It may be possible to adapt this proof to prove the above theorem for the setting in which a set 
is not allowed to be included more than once. This, however, appears to be nontrivial and instead 
we next give a proof, in Section ISTTl that uses the LP to reduce the problem to a regular set cover 
problem with a modified set system whose primal and dual VC dimensions are at most 0{S) and 
0((5*), respectively. 

3.1 Multi-cover without repetition 

Geometric intuition. Imagine we have a set of points and a set of disks = {ri, . . . ,rm} (i.e., 
the ranges) in the plane. We solve the LP for this system. This results in weight assigned to each 
disk, such that the total weight of the disks covering a point p G P exceeds (or meets) its demand 
d(p). We add another dimension (we are now in three dimensions), and for each i = l,...,m 
translate the disk rj G 9^ to the plane z = i. Let 9"' denote the resulting set of m two-dimensional 
disks that "live" in three dimensions. Observe that the projection of 9"' to the xy plane is 9". Every 
point pj € P is now a vertical line tj (parallel to the z-axis), and we are asking for a subset X of 
9"', such that every line ij stabs at least d(pj) disks of X. The fractional solution for the original 
problem induces a fractional solution to the new problem. The next step, is to break every line ij 
into segments, such that the total weight of the disks of 7' intersecting a vertical segment is at least 
1 (and at most 2). Let L' be this resulting set of segments. Consider the "set system" S = (L', ?"'), 
and its associated set cover instance of the disks of "J' so that they intersect all the segments of L' . 
It is easy to verify that any solution of this set cover problem, is in fact a solution to the original 
multi-cover problem, and vice versa (up to small constant multiplicative error, say 2). We know how 
to solve such set-cover problems using standard tools. The key observation is that the projection 
of (L',?'') on to the plane yields the original range space. Similarly, projecting {L',3^') on to the 
z-axis results in a range space where the points are on the real line and the ranges are intervals. 
Since the range space {L',3^') is the intersection of two range spaces of low VC dimension, it has 
low VC dimension. This implies that the set-cover problem on {L',S'') has a good approximation 
[BG95| and this leads to a good approximation to the original multi-cover problem on S. 
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More formal solution. Consider a fractional solution x to the LP associated with Z. If any set 
G 9" satisfies rcj > 1/4 then we add rj to our solution. There can be at most 4^jXj = 4f such 
sets, so including them does not harm our goal of a solution with 0(f) sets. We now work with the 
residual instance and hence we can assume that the fractional solution has no set with > 1/4. 

Now, assume that we have fixed the numbering of the ranges of 9" = {ri, . . . , rm}, and consider 
the fractional solution, with the value Xj associated with r^, see Eq. In particular, for a point 
p G P, consider the linear inequality 

^ Xi> d(p) . 

This inequality holds for the fractional solution. We split this inequality into 0(d(p)) inequalities 
having 1/2 on the right hand side. To this end, scan this inequality from left to right, and collect 
enough terms on the left-hand side, such that their sum (in the fractional solution) is larger than 
1/2. We will write down the resulting inequality, and continue in this fashion until all the terms of 
this inequality are exhausted. 



Formally, let Uq = U 



p G r.j I be the sequence of indices of the ranges participating in the 
above summation, where U and Uq are sorted in increasing order. For £ > 1, let be the shortest 
prefix of C/^-i such that X^igv^ ^« — -'-Z^' ^^'^ largest number (i.e., index) in V^, and 

let U£ = {Ue^i \ Vi). Since each Xj < 1/4 we have that J2ieVe ^« ^ -'-/^ -|- 1/4 < 3/4. We stop when 
"^ieUe < 1/2 for the first time. This process creates some h inequalities of the form 

^x, > 1/2, 

i£Ve 

for £ = 1, . . . , h. We have h > d(p) inequalities from the fact that X^j-pgr^ ^ <^(p) ^iid by our 
observation that X^j^vi Xi < S/A. 

We next describe a new set system (P',9"), derived from this construction of inequalities, such 
that a set cover solution to the new system implies a multi-cover solution to the original system, 
and the new system has small VC dimension. 

The new set system (P', 3") is defined as follows. For each point p which was processed as above, 
we create h copies of it, one for each V^. Each such copy of p corresponds to an interval I = [a, 13], 
where a is min^gy^ i, and (3 is maxjgv^ i. So p has h such intervals associated with it, say Ii, . . . ,Ih- 
We generate h new pairs from p, namely, Q{p) = {(p, Ii) , . . . , (p, Ih)}- 

We set P' = UpQ(p), and 7 = jrj G where 

(p,/)gP' pGrjandiG/}. (2) 



ri 

Note that \fi\ = |rj|, and it can be interpreted as deciding, for each point p G which one of its 
copies should be included in f^. 

The following two claims follow easily from the construction. 

Claim 3.2 For the set cover instance defined by (P',3") there is a fractional solution of value 
2^.x, <2f. 

Claim 3.3 An integral solution of value (3 to the set cover instance (P',3") implies a multi-cover to 
the original instance of cardinality at most (3. 

We need the following easy lemma on the dimension of intersection of two range spaces with 
bounded VC dimension. 
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Lemma 3.4 ( |Har08| ) Let S = {X,TZ) and T = {X,TZ') be two range spaces of MC- dimension 6 
and 5' , respectively, where 5,5' > 1. Let 7^ = |r PI r' r € 7^, r' € 7^' | . Then, for the range space 
S= {X,ii), we have that 5$) = 0{5 + 5'). 

Observation 3.5 IfS = {X,TZ) hasVC dimension 5, andM. C TZ, then theVC dimension o/(X, M) 
is bounded by 5. 

The crucial lemma is the following. 

Lemma 3.6 The VC dimension of the set system (P',9") is 0(5) and the VC dimension of its dual 
set system is 0{5*). 



Proof: We define two set systems ( P',9"j and (P',9") as follows. 3' 



where 



and 9" 



f; = {(p,/)GP' 

G jj , where fj = | (p, /) G P' z G / 1 . 



P G Ti 



Note that fi = ri D ri (see Eq. Therefore (P',9") is formed by the intersection of ranges 

(P', 9") with ranges of (P', 3"). Therefore the VC dimension of (P', 9") is bounded by 0^5 + 5^ where 

5 and 5 are the VC dimensions of (P', 3') and (P', 9") respectively, by Lemma [3^ and Observation ESI 
We observe that the set system (P',9") has the same VC dimension as that of (P,3") since we only 
duplicate points. The set system (P',3") has constant VC dimension 5 = 3 since it is the intersection 
system of points on the line with intervals. 

The second part of the claim follows by a similar argument. Consider the dual range spaces of 
^P', 3"^ , (P', 9") , and ^P', 9"^ , respectively. The ground set of these range spaces can be made to be 
H. We have the following: 



• X* = l^y, Mj, the range space dual to yP',3'j, has for any point (p,/) G P' a range that 

contains all the ri that contains p. It is therefore just the dual range space ioX = (P, 7), 
and it has VC dimension 5* . 

• X* = M), the range space dual to (P',9"), for every (p, /) G P', has the range containing 
all the sets such that i G /. As such, X* has a constant VC dimension. 

• X* = , the range space dual to ^P', 3^ , for every (p, /) G P', has the range containing 
all the sets r, such that i G / and p G r^. 

We have that X* is the range space contained in the intersection of range spaces X* and X*. 
Lemma [3^ and Observation 13.51 imply that the VC dimension of X* is 0{5*). ■ 

Now we apply the known results on the integrality gap of the LP for set cover as discussed in 
Section [2?3l These results imply that for the set system {P',3^) there is an integral set cover of 
value 0((5*flogf) (here we use Claim [321 and Lemma [321) ■ Prom Claim [331 there is a multi-cover 
for the original instance of the desired size. This completes the proof of the theorem. We observe 
that the algorithm is in fact quite simple. After solving the LP, pick each range rj independently 
with probability min{l, cxj} where c = a ■ 5* logf for a sufficiently large constant a. With constant 
probability this yields a multi-cover. 
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3.2 Multi-cover in spaces with bounded VC dimension when allowing repetition 



We consider the case where sets in 3^ are allowed to be picked multiple times to cover a point. For 
this purpose we use relative approximations. The standard definition of relative approximation is 
the dual of what we give below. 

Let a, (/) > be two constants. For a set system X = (P,^), recall from Definition 12.11 
that #(pn9") denotes the number of sets in 3^ that contain the point p. A relative {a^cj))- 
approximation is a subset X C 3" that satisfies 

for each p € P with #(p n 9") > (j) • It is known [LLSOlj that there exist subsets with this 
cS 1 

property of size log — , where c is an absolute constant, and 5 is the VC dimension of the dual 

set system of (P,?'). Indeed, any random sample of that many sets from 9" is a relative {a,(p)- 
approximation with constant probability. To guarantee success with probability at least 1 — q, one 

needs to sample I (5 log — + log - ) elements of X for a sufficiently large constant c ILLSOII . 
a^<p \ cp qj 

To apply relative approximation for our purposes we let be a large integer such that Nxi is 
an integer for each range rj (since the Xi are rational such an N exists). We create a new set system 
{P,3^') where 5"' is obtained from 3^ by duplicating each range G Nxi times. Thus I^T'I = Nf. 
From the feasibility of x for the LP we have that #(p n 3'') > iVd(p) > A^/d(p) // for each p G P. 

Now we apply the relative approximation result to (P,3"') with cp = 1/f ^ — 1/2 to obtain 
a set X C 9"' such that \X\ = 6(5*f logf) and with the property that for each p S P, 

#(pn J') ^ #(pnx) 



2\y\ - \X\ 
We have 

«p n X) > m . mnn ^ m , ^ . „p) ^ , .(p, , 

as desired. 

Note that X is picked from H' which has duplicate copies of sets from 9". Recall that the 
algorithm, from the previous section (which is for the variant without repetition), picks each range 
ri independently with probability min{l, cxi ■ 6* log f }; and this yields a feasible multi-cover without 
repetitions. It may be possible to analyze this algorithm (i.e., without repetitions) directly by a 
careful walkthrough of the proof for relative approximations. 



4 Multi-cover for Halfspaces in 3d and Generalizations 

In this section, we show that improved approximations can be obtained for specific classes of set 
systems induced by geometric shapes of low complexity. In particular, we describe an 0(1) ap- 
proximation for the multi-cover problem when the points are in H^, and the ranges are induced 
by halfspaces. The main idea, of using cuttings, extends also to other nice shapes. We outline the 
extensions and some applications in Section [5l 
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4.1 Total demand, Sampling, and Residual demand 

We develop some basic ingredients that are useful in randomly rounding the LP solution. These 
ingredients apply to a generic multi-cover instance, not necessarily a geometric one, however we use 
the notation of points and ranges for continuity. 



Lemma 4.1 Given a multi-cover instance T 
more than the total demand dj(P). 



[P,3'), one can compute a cover for 2 of size no 



Proof: Indeed, scan the unsatisfied points of P one by one, and for each such point p, add to the 
solution d(p) ranges that cover it, picked arbitrarily. Clearly, the ranges that are picked satisfy all 
the demands, and the number of ranges picked is at most X]pd(p) = dj(P). ■ 

Given an instance of multi-cover I = (P,9^) and a feasible fractional solution x, a cx-sample 
for a scalar c is a random sample of J", formed by independently picking each of the ranges Vi £ 3' 
with probability min{l,ca;j}, where Xi is the value assigned to rj by the fractional solution. (For 
the i with cxj > 1, so that i is chosen with probability one, we will simply assume that such choices 
have been made, and the demand removed; that is, we assume that hereafter that Xi < 1/c. Since 
the number of such i is at most cf, this step does not affect our goal of obtaining an output cover 
with 0(f) sets.) 

Lemma 4.2 Let c> 4 be a constant and let I = (P, 5") be a multi-cover instance with an LP solution 
satisfying Xi < 1/c for all i. If % is a cx-sample and p G P zs a point with demand d = d(p), then 



Pr 



p is not fully covered by 3i 



Pr 



#(pn3i) < d 



< exp 



and E 



dres(p,3i) 



< exp(--d 



Proof: Let Xi be the indicator variable which is equal to one if the cx-sample includes the range 
Ti G If, and is zero otherwise. Let Y = #(p n 3i) = Yli-p^n observe that ^ = 'E\Y] > cd using 
the facts that x is a feasible solution to LP, and xi < 1/c for all i. For j G [0,d], we apply the 
Chernoff inequality |MR95j and use the fact that c > 4 to obtain: 



Pr 



#(pn3?) <d-j 



< Pr 



y < Ml - (c - i)/c - jV/i) 

H 3 



/ fifc-lj 



< exp(-^--jl <exp( -^d-^j 



The first statement of the lemma follows by substituting j = 1 and observing that the desired bound 
follows, and the second follows by using the fact that, for a random variable Z taking non-negative 
integral values, that 'E[Z] = ^^^QPr[Z > k]. This implies 



E 



dres(p,3i) 



l<i<d 



exp(--d 



#(p n 3i) < d - j 

3 



E 

l<i<d 



exp 



l<j<d 

c . 



"4^ - 4^ 



-jj <exp^ exp(3/4)-l 



< exp( --d 



as claimed. 
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In the following, for t > 1, let 



Pt= {pe P I t<d(p) <2t }. 



The lemma below implies that if the number of points in the set system is "small" then the 
multi-cover problem can almost be solved in one round of sampling. 

Lemma 4.3 Suppose there is a probability distribution on a collection of multi-cover instances such 
that an instance I = (P, 9") chosen from the distribution satisfies, for any t > 1, that 



E 



<V -K^ 



where K and V are fixed parameters of the distribution. Then there is a value c depending on 
K , such that for any feasible fractional solution x to I, a cx-sample 01 results in expected total 
residual demand drcs(P)3i) < V; here the expectation is with respect to the randomness of 2 and the 
independent randomness of the cx-sample. 

Proof: Let 01 he a cx-sample of 3^ for fixed c > 4 + 4 log K. Let X be the subset of S' with all 
ranges having Xi > 1/c. Since 3i \ X is also a cx-sample of 1\ X, we assume hereafter that X is 
empty; the result for general X follows by application of the result to T\ X. 

By applying Lemma [42l to the induced range space (P^,?"), we have 



E 



drcs(Pt; 3?) 



< Ei 



j;E3,[dre.(p,3?)] 

pGPt 



< Ej 



\Pt\exp(^-^t 



Ej 



exp(--i 



< VlCexp(--t) < Vexp{-t{c/4:-\ogK)) < V exp{-t) 



Then, by linearity of expectation, we have 



E 



dres(P, Oi) 



EE 

i=0 



dres(P2M ^) 



< ^yexp(-2*) < V. 



i=0 



Thus, after cx-sampling, the residual instance has total expected demand bounded by V, as claimed. 



4.2 Clustering the given instance 

The key observation to solve the multi-cover problem in our settings is Lemma [4?3l as it provides 
a sufficient condition for an 0(1) approximation. Of course, it might not be true (even in low 
dimensional geometric settings) that the number of points (i.e., the total residual demand) is small 
enough, as required to apply this lemma. We preprocess the given instance via an initial sampling 
step and then employ a clustering scheme that partitions the points into regions; we argue that 
these regions and an induced multi-cover instance on them satisfies the conditions of the lemma. 

The depth of a simplex A in a set of weighted halfspaces is the minimum depth of any point 
inside A, see Definition 12.11 

To perform the aforementioned clustering, we will use the shallow cutting lemma of Matousek 
[Mat92j . We next state it in the form needed for our application, which is a special case of Theo- 
rem ETH 
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Lemma 4.4 Given a set H of weighted halfspaces in IR^, with total weight W , there is a randomized 
polynomial-time algorithm that generates a set T of simplices, called a (l/4Ty)-cutting, with the 
following properties: the union of the simplices covers H^; the total weight of the boundary planes 
of 3^ intersecting any simplex ofT is bounded by 1/4; and finally, for any t >0, the expected total 
number of simplices of depth at most t is O(VFt^). (Here the expectation is with respect to the 
randomness of the algorithm.) 

4.2.1 The algorithm 

Given an instance of multi-cover 2 = (P,9") of points and halfspaces in IR'^, our algorithm first 
computes the fractional solution to the LP induced by X, yielding weights Xi. Next, for (3 an 
absolute constant in (0, 1/4) to be specified later, we put in the set X all the ranges rj with Xi > (3. 
Let (P', ff'') = (P, \ X. Let f = 'Yliri&'3\x total weight of the remaining ranges. 

The remainder of the algorithm uses a auxiliary abstract multi-cover instance derived using 
cuttings, as described next. 

Using the weights Xj, we build a (l/4f')-cutting F for 3"'. This induces an abstract multi-cover 
instance (T,S''), where a simplex A € F is covered by halfspace h £ 3'' only if the interior of A is 
contained inside h and it does not meet the boundary plane of h. The demand cl(A) is defined to 
be maXpgpp^ cli-es 

A feasible solution to (F, 'J') is also, by construction, a feasible solution for the original instance 
I. Furthermore, any feasible fractional solution for I can be transformed into a feasible fractional 
solution for (F, 9^'), at the cost of a constant factor. Indeed, the weights Xi give a feasible fractional 
solution to I\X, and so the depth of A is at least d(A) — 1/4, where A "loses" at most weight 1/4 
of depth due to halfspaces whose boundary planes cut A. It follows that if the depth is measured 
with respect to weights Xi = 2xi, the new depth of A (i.e., the point with minimum cover in A) 
is at least 2d(A) — 1/2 > d(A). That is, the weights Xi give a feasible fractional solution to the 
multi-cover instance (F,3''). Note that since /? < 1/4, the weights Xi satisfy Xi < 1, for all i. 

The remainder of the algorithm is to apply the approach implied by Lemma 14.31 we find a cx- 
sample Ji, with c to be determined; this induces a residual multi-cover problem (F, 5"') which we 
solve using the simple technique of Lemma 14.11 Letting U denote the resulting combined solution 
to (F,3"'), we return C/ U X as a cover for the original multi-cover problem. 

The analysis of this algorithm is the proof of the following result. 

Theorem 4.5 Let 2 = (P, 3^) be an instance of multi-cover formed by a set P of points in IR'^, and 
a set 9" of halfspaces. Then, one can compute, in randomized polynomial time, a subset of halfspaces 
of 3' that meets all the required demands, and is of expected size 0(f), where f is the value of an 
optimal fractional solution to LP. 

Proof: We described the algorithm above, except for the values of c and 13. 

By Lemma 14.41 the expected number of simplices in the cutting F of demand at most t is 
0{Wt'^^, where W = f < f(X), which implies that Lemma 143) can be applied, with V = f(X), K an 
absolute constant, and using the weights Xi. Since Xi < 2f(Z), the expected size of U is at most 
(c-|- 2)f (X), using the absolute constant value of c used in this application of Lemma [43l Observing 
that \X\ < f(X) //3, and taking P = l/2c to allow the cj;-sample probabilities cxi to be less than 
one, we have that the returned solution [/ U X to X has expected cardinality at most (3c -|- 2)f(X), 
which is 0(f(X)). 

The only non-trivial step in terms of verifying the running time is for computing the cutting 
and Lemma Hm guarantees the running time. ■ 
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Remark 4.6 The shallow-cutting lemma (Lemma I4.4|) is shown via a random sampling argument, 
and our rounding algorithm is also based on random sampling, given the cutting as a black-box. 
One could do a direct analysis of random sampling by unfolding the proof of the cutting lemma. 
However, the indirect approach is easier to see and highlights the intuition behind the proof. 

5 Generalizations and Applications 

We now examine to what extent the result derived for covering points in by halfspaces generalizes 
to other shapes. 

5.1 Well behaved shapes 

We are interested in set systems (P,9") where 9" is a set of "well-behaved" shapes such as disks or 
fat triangles. As we remarked already, it is shown in |CV07| that the existence of good e-nets for 
such shapes can be derived from bounds on their union complexity. For example, it is shown that if 
y is a set of fat triangles in the plane then there is an O(loglogf) approximation for the set cover 
problem. For fat wedges one obtains an 0(1) approximation. Here we show that union complexity 
bounds can be used to derive approximation ratios for the multi-cover problem that are similar to 
those derived in |CV07j for the set cover problem. Following the scheme for halfspaces, the key tool 
is the existence of shallow cuttings. To this end we describe some general conditions for the shapes 
of interest and then state a shallow cutting lemma. 

Let y be a set of n shapes in K'^, such that their union complexity for any subset of size r 
is (at most) U{r), for some function U{r) > r. Similarly, let 0[r'^) be the upper bound on total 
complexity of an arrangement of r such shapes. 

Let X be a subset of K'^. We assume that given a subset S ^ 3^, one can perform a decomposition 
the faces of the arrangement ^(S) that intersects X into cells of constant descriptive complexity 
(e.g., vertical trapezoids), and the complexity of this decomposition is proportional to the number 
of vertices of the faces of ^(S) that intersects X. Finally, we assume that the intersection of d 
shapes of 3^ generates a constant number of vertices. 

One can then derive the following version of Matousek's shallow cutting lemma. We emphasize 
that this lemma is a straightforward (if slightly messy) adaption of the result of Matousek. A proof 
is sketched in Appendix lAl 

Theorem 5.1 Given a set If of "well-behaved" shapes in IR'^ with total weight n, and parameters 
r and k, one can compute a decomposition of space into 0{r'^) cells of constant descriptive com- 
plexity, such that total weight of boundaries of shapes of 3' intersecting a single cell is at most n/r. 
Furthermore, the expected total number of cells containing points of depth smaller than k is 



where U{t) is the worst-case combinatorial complexity of the boundary of the union of I shapes of 



Using the same scheme as that for halfspaces we can derive approximation ratios for the multi- 
cover problem for shapes that have the property that lA{n) is near-linear in n. An approximation 
ratio of 0(^/(opt) /opt) easily follows, but in fact, by using the oversampling idea of Aronov et al. 
[AESOQj . we can improve this to 0(log(Z//(opt) /opt)). We use the shallow cutting lemma as a black 
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box, and hence our argument is arguably slightly simpler than then one in |AES09| and our result 
can be interpreted as a generalization. 



Theorem 5.2 Let 2 = (P, 9") be an instance of multi-cover formed by a set P of points in TR!^, and 
a set 3^ of ranges. Furthermore, the union complexity of any i such ranges is (at most) U{£), for 
some function U{£) > i. Then, one can compute, in randomized polynomial time, a subset of ranges 

of 9' that meets all the required demands, and is of expected size O^flog^^^, where f is the value 
of an optimal fractional solution to LP. 



Proof: As before, we compute the LP relaxation, and take all the ranges that the value of Xi > /3, 
where jS = a/ log for some sufficiently small constant a. Next, we compute a (l/4f)-cutting T 
of residual system (P',9"'). Using Theorem 15.11 with parameters r = 4f , n = f and k = t, there are 
at most 

cells, with depth at most t. In particular, this bounds the number of cells in the cutting with depth 
in the range t — 1 to t. We pick a random sample 3i of (expected) size h = O^f log from 9"', 

by performing a CT-sample from 9"', where c = O^log • Arguing as in Lemma [4.21 the expected 

residual demand for a cell of T with demand t is t exp(— ct/4). Therefore, the expected total residual 
demand in (F, 9"') \ 3i is 



off;-p(4t) (t+i)^+^f^j=o(f) 



Using Lemma [411 the residual multi-cover instance (r,9"') \ 3i has a cover of expected size 0(f). 
Thus, we have shown that the original multi-cover instance has a cover of expected size 0(f//? + 



/i + f) = 0(^flog^j 

Applications: The above general result can be combined with known bounds on U{n) to give 
several new results. We follow |CV07l lAESQQj who gave approximation ratios for the set cover 
problem using a similar general framework; we give essentially similar bounds for the multi-cover 
problem. All the instances below involve shapes in the Euclidean plane. 

• 0(1) approximation for pseudo-disks, fat triangles of similar size, and fat wedges. 

• 0(log log log f ) approximation for fat triangles (which also implies similar bounds for fat convex 
polygonal shapes of constant description complexity). 



0(logQ(f)) approximation for regions each of which is defined by the intersection of the non- 
negative y halfplane with a Jordan region such that each pair of bounding Jordan curves 
intersecting at most three times (not counting the intersections on the x axis). Here a{n) is 
the inverse Ackerman function. 



5.2 Unit Cubes in 3d 

We also get a similar result for the case of axis-parallel unit cubes. 

In [CV07| an 0(1) approximation is also shown that for the problem of covering points by unit 
sized axis parallel cubes in three dimensions. There is a technical difficulty for this case. Although 
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it is known from [BSTY98j that the combinatorial complexity of the union of n cubes is 0{n), the 
same bound is not known for the canonical decomposition of the exterior of the union as required by 
our framework. The same difficulty is present in [CV07| and they overcome this by taking advantage 
of the fact that all cubes are unit sized. The basic idea is to use a grid shifting trick to decompose 
the given instance into independent instances such that each instances has cubes that contain a 
common intersection point. For this special case one can show that the canonical decomposition of 
the exterior of the shapes has linear complexity. This suffices for the framework in |CV07j . For our 
framework we need a cutting. 

Lemma 5.3 Let S he a set of n axis-parallel unit cubes in three dimensions, all of them containing 
(say) the origin. Then, one can decompose the arrangement of A{S) into a canonical decomposition 
of axis parallel boxes, such that the complexity of decomposing every face is proportional to the 
number of vertices on its boundary. 

Proof: First we break the arrangement into eight octants by the three axis planes {xy, yz and 
xz planes). We will describe how to decompose the arrangement in the positive octant, and by 
symmetry the construction would apply to the whole arrangement. 

So, let / be a 3d face of the arrangement (when clipped to the positive octant). Let I be the 
cubes of S that contain /, and similarly, let B be the set of cubes of S that contribute to the 
boundary of /, but do not include / in their interior. As such, we have that 



(If the set I is empty, we will add a fake huge cube to ensure / is bounded.) Now, the first term 
is just an axis-parallel box. Intuitively, the second term (the "floor" of /) is a (somewhat bizarre) 
collection of "stairs". Note, that any vertical line that intersects /, intersects it in an interval. In 
particular, let g the top face (in the z direction) of /, and observe that, since all the cubes of S 
contain the origin, it must be that any line that intersect / must also intersect g. As such, let us 
project all the edges and vertices of / upward till the hit g. This results in a collection W of (interior) 
disjoint segments that partition (the rectangular polygon) /. We perform a vertical decomposition 
of the planar arrangement formed by yi(VF) (including the outer face of this arrangement, which is 
g). This results in 0(|/|) collection of (interior) disjoint rectangles that cover g, where |/| is the 
number of vertices on the boundary of /. Furthermore, for such a rectangle r, there is no edge 
or vertex of /, such that their vertical projection lies in the interior of r. Namely, we can erect a 
vertical prism for each face of the vertical decomposition olAiW), till the prism hits the bottom 
boundary of /. This result in a decomposition of / into 0(|/|) disjoint boxes, as required. ■ 

Lemma [531 implies that an the arrangement A{S), can be decomposed into (canonical) boxes, 
in such a way that the number of boxes of certain depth t, is proportional to the number of vertices 
of A{S) of this depth. This implies that we can apply the shallow cutting lemma to S (we remind 
the reader that all the axis-parallel unit cubes of S contain the origin). 

This is sufficient to imply 0(1) approximation to multi-cover. Indeed, let X = (P,?") be the 
given instance of multi-cover, where is a set of unit-cubes in three dimensions. Let G be the unit 
grid, and for any point q € G, let be the set of cubes of 9" that contains p (for the simplicity 
of exposition, we assume that every cube of 3" is contained in exactly one such set, as this can be 
easily guaranteed by shifting G slightly). Next, solve the LP associated with X, and associate a 
point p € P with q G G, if the depth of p in 9"q is at least 1/8 (if p can be associated with several 
such instances, we pick the one that provides maximum coverage for p). Let Pq be the resulting 
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set of points. Thus, for any point in q € G, there is an associated instance of muhi-cover (Pq,9"q). 
Clearly, a constant factor approximation for each of these instances, would lead to a constant factor 
approximation for the whole problem. 

Now, IFq is made of cubes all containing a common point, and as such Lemma [531 implies that 
shallow cutting would work for it. In particular, we can now apply the algorithm of Theorem 14.51 to 
this instance, and get a constant factor approximation (here, implicitly, we also used the fact that 
the union complexity of n axis-parallel unit cubes is linear). This implies the following theorem. 

Theorem 5.4 Let 2 = (P, 9") be an instance of multi-cover formed by a set P of points in M^, and 
a set "J of axis-parallel unit cubes. Then, one can compute, in randomized polynomial time, a subset 
of cubes of 3^ that meets all the required demands, and is of expected size 0(f), where f is the value 
of an optimal fractional solution to LP. 

6 Conclusions 

We presented improved approximation algorithms for set multi-cover in geometric settings. Our 
key insight was to produce a "small" instance of the problem by clustering the given instance. This 
in turn was done by using a variant of shallow cuttings. We believe that this approach might be 
useful for other problems in geometric settings. 

An interesting open problem, is to obtain improved algorithms for the set cover and the set 
multi-cover problems in geometric settings when the sets/shapes have costs associated with them 
and the goal is to find a cover of lowest cost. Can the results from [Cla93l IBG951 ICV07| and this 
paper be extended to this more general setting? 

Recently, Mustafa and Ray [MR09| gave a PTAS for the problem of covering points by disks in 
the plane; their algorithm is based on local search. It would be interesting to see if this algorithm 
can be adapted to the multi-cover problem. 
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A A shallow cutting lemma for "nice" shapes 

In this section, we prove Theorem 15.11 a variant of the shallow cutting lemma of Matousek in a 
slightly different setting. We include the details for the sake of completeness, which are not hard in 
light of Matousek's work |Mat92| . Our description is somewhat informal, for simplicity. The family 
of shapes that we consider needs to satisfy the assumptions outlined in Section O 

Building (l/r)-cuttings. When computing cuttings, one first picks a random sample of size 
r of the objects of and computes the decomposition Ai^\{^) of the arrangement of the random 
sample. For a cell A in this decomposition, let cl(A) be the list of shapes of 3' whose boundaries 
intersect the interior of A. If |cl(A)| < n/r then it is acceptable, and we add it to the resulting 
cutting. 

Otherwise, we need to do a local patching up, by partitioning each such cell further. Specifically, 
let tA = rcl(A) /(n/r)] be the excess of A. We take a random sample Ji^ of size 0{t/\log{t/\)) 
from cl(A). With constant probability, this is a l/t^^-net of cl(A) (for ranges formed by our 
decomposition). We verify that it is such a net, and if not, we resample, and repeatedly do so 
until we obtain a l/t/\-net. To do the verification, we build the arrangement of D^a inside A, and 
compute its decomposition, and check that all the cells in this decomposition intersect at most n/r 
boundaries of the shapes of 1. Let dcmp(A) denote this decomposition of A (if A has excess at 
most 1, then we just take dcmp(A) to be {A}). Clearly, the set 

IJ dcmp(A) 

forms a decomposition of IR*^ into regions of constant complexity, and each region intersects at most 
n/r boundaries of the shapes of 3^. 

It is well known that the complexity of the resulting cutting is (in expectation) 0{r'^) |CF90| 

Let S denote the resulting cutting. 



Size of cutting at a certain depth. Here we are interested in the number of cells in the 
arrangement >l||(3i) that cover "shallow" portions of ^1(9"). Formally, the depth of a point p G IR"^, 
is the number of shapes of 3^ that cover it. Let f<k{n) denote the maximum number of vertices of 
depth at most k in an arrangement of n shapes. Clarkson and Shor |CS89j showed that f<k{n) = 
0(^k'^U{n/k)) . Specifically, we are interested in the number of cells of C that contain points of depth 
at most k. The kth level is the closure of all the points on the boundary of the shapes that are 
contained inside k shapes. 

Now, the expected number of vertices of A{^) that are of depth at most k in A{3') is 



o(^(^-^y k'^U{n/k)^ =oU^^ U{n/k) 



since for a given vertex of A{3^) of depth at most fc, the probability that all d shapes that define it 
will picked to be in % is 0((r/?7,)'^). This unfortunately does not bound the number of cells in the 
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decomposition of that contain points of depth at most k, since we might have cells that cross 

the kth level. 

So, let X C TR'^ be a fixed subset of space, and let be the number of cells of that 

intersect X. Let x{r) denote the maximum value of over all samples of size r. Similarly, 

let x*'{Jl) denote the number of cells in that intersect X and have excess more than t (i.e., 

there are at least t ■ n/r shapes intersecting this cell). 

Chazelle and Friedman [CF90| showed an exponential decay lemma stating that e[x*(3?)] = 
0(2"* B[x{^)]) . We comment that, in fact, one can prove directly from the Clarkson-Show technique 
a polynomial decay lemma, which is sufficient to prove the shallow-cutting lemma. This polynomial 
decay lemma is implicit in the work of de Berg and Schwarzkopf |dBS95j although it was not 
stated explicitly (it also made a stealthy appearance in Clarkson and Varadarajan work |CV07| . 
but [dBS95| seems to be the earliest reference). 

Lemma A.l (Polynomial decay lemma.) For t > 1, let % he a random sample of size r from 
"J, and let c>l be an arbitrary constant. Then E[x*(DI)] = 0{x{r) /f^) . 



Proof: By the Clarkson-Shor technique |CS89l lCla88j . we have that 



E 



Aeyi||(K) 



O 



^yn<n)=o[(^yx{r)). 



In particular, if there are x^{'X) cells in ^|(3?) with conflict-list of size larger than t{n/r), then they 
contribute to the left size of the above equation the quantity x^{Ol){t{n/r)Y. We conclude that 

B[x\mHr)r]=0(l^yx{r)), 
which implies that E[a;*(3i)] = 0{x{r)/t'^), as claimed. ■ 

Lemma A. 2 The expected number of cells in the {l/r)-cutting Q of 3^ that contain points of depth 
at most k is bounded by 




Proof: If a cell A of ^|(3?) has excess t, and it intersects the kth level, then all its points 
have depth at most k + t(n/r). The expected number of vertices of ^|(3i) of depth at most 
a{t) = k + t{n/r) is 



n 



n J \Ci{t) J j 

which also (asymptotically) bounds the number of cells in ^|(3?) having depth smaller than a{t). 
Let Xt denote the number of cells with excess t (or more) with depth at most a{t). Setting c = 0{d), 
we have by the polynomial decay lemma, that 



t^n J \a{t) J I 

Now, the number of cells of the cutting S that have points with depth at most k is bounded by 

/ oo \ 



y = o[Y,Xf{tiogtr 
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Thus, we have 

by setting c to be sufficiently large. ■ 
The above proves Theorem 15.11 by using replication to represent weights. 
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